Decay of spin-1/2 field around Reissner-Nordstrom black hole 
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To find what influence the charge of the black hole Q will bring to the evolution of the quasinormal 
modes, we calculate the quasinormal frequencies of the neutrino field (charge e = 0) perturbations 
and those of the massless Dirac field (e 7^ 0) perturbations in the RN metric. The influences of Q, e, 
the momentum quantum number I, and the mode number n are discussed. Among the conclusions, 
the most important one is that, at the stage of quasinormal ringing, the larger when the black hole 
and the field have the same kind of charge (eQ > 0), the quasinormal modes of the massless charged 
Dirac field decay faster than those of the neutral ones, and when eQ < 0, the massless charged Dirac 
field decays slower. 
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It is well known that black holes and neutron stars have their characteristic "sound": quasinormal modes, which 
\Q ■ were originally found to play a prominent role in the gravitational radiation theoretically and were hoped to be 
detected in the gravitational experiments (see 1, 2] and references therein for review). They have been extensively 
studied recently for their possible connection with the quantum gravity through an observation made by Dreyer in 
Ref. j3j, where the author fixed the free parameter 7 of Loop Quantum Gravity(LQG) and suggested the gauge group 
of LQG may be SO(3) rather than SU{2). This is remarkable^ and surely makes a good example of utilizing the 
' macroscopic behavior of black holes to deduce their microscopic nature. 

Besides the application in LQG, the quasinormal modes are also used in many other fields, especially the Anti-de 
Q\ • sitter /Conformal Field Theory (AdS/CFT) correspondence^. In Ref. 6], from the AdS/CFT correspondence, the 
quasinormal modes of the AdS black hole are used to determine the relaxation time of a field perturbation, i.e., the 
C*") 1 imaginary part of the quasinormal frequency is proportional to the inverse of the damping time of a given mode. 

Hod and Piran0 considered the behavior of a charged scalar field perturbation and found that, its quasinormal 
modes will dominate the radiation in the late time evolution for their slower decay than the neutral ones. Konoplya, 
in his workp|, first investigated a complex (charged) scalar field perturbation through calculating its quasinormal 
$_i , modes in the Reissner-Nordstrom (RN), RNAdS and dialton black holes and found that, on the contrary, the neutral 
perturbations will decay slower than the charged ones at the stage of quasinormal ringing. In this paper, we turn 
I to treat the spin 1/2 Dirac particles, in the presence of the RN black holes with charge Q. Calculating quasinormal 
frequencies of the neutrino field (charge e = 0) perturbations and those of the massless Dirac field (charge e ^ 0) 
perturbations in the RN metric, we found that when the product eQ is positive, the quasinormal modes of the massless 
Dirac field decay faster than those of the neutral ones and when eQ is negative, the perturbations of the massless 
Dirac field decay slower than those of the neutral field. The detailed discussions can be found in the following sections. 

In Sec. [n] in the RN metric, we will calculate the quasinormal frequencies of the neutrino field perturbations and 
those of the massless Dirac field perturbations. In Sec. II I II with some figures and tables, we show clearly what 
influence the charge of the black hole brings to the behavior of the quasinormal modes, including the energy and the 
decay rate. Sec. IIVI is a summary with some suggestions about the future research. 

II. QUASINORMAL MODES IN RN METRIC 

We shall calculate the quasinormal frequencies of the massless Dirac field perturbations and the neutrino field 
perturbations in the RN metric 

ds 2 = -(l-^E + ^dt 2 +(l-'^ + ^j 1 dr 2 +r 2 (d9 2 + sin 2 9dcf> 2 ) . (1) 
The wave equation can be written as 

h a e c r(d fl + r fi + e A^ = o. (2) 
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In Eq. (J2J, 7° are the Dirac matrices with the forms 



where tr* are the Pauli matrices. e a M is the inverse of the tetrad e„° defined by the spacetime metric 

9nv = VabefeJ 1 , t] ah = diag(-l, 1, 1, 1). 

are the spin connection given by 

where r" are the Christoffel symbols. And, e^4 M is the electromagnetic potential of the black hole. 

In the following discussions, we adopt the natural unit, h = c = G = 4tteq = 1. It means that all physical quantities, 
such as mass, length, charge, time, energy, and so on, are dimensionless. The following replacing relations can be 
used to restore all of them: 

/ ch\ m [ch 

Mass : m — ► — m = — , m„ = \ — , 
\ G J m p V G 

HG\~ 1/2 , I , fhG 



Length : / — I -g-J I = — , Z p — 
lime : i — m — — I t = — , i p = 




Cm 

e 



Charge : e — > (47reocft) e = — , e p = \J A-keqcH, 
c 5 h\~ 1/2 „ £ „ /c^fi 



Energy : £ - ( — j £ - — , E p y f . 



For the RN spacetime, the electromagnetic potential eA^ can be written as 



eA M = 1^,0,0,0), (>!) 



r 



where Q is the charge of the RN black hole, and e is zero for neutrinos and nonzero for Dirac particles. In the following 
we shall see that we can just consider the sign of the product eQ without taking into consideration the sign of each 
quantity respectively. Therefore, for convenience we shall take Q to be always positive and only switch the sign of e, 
when we are studying the influence of the charge of the black hole on the quasinormal frequencies of the Dirac field 
perturbations. 

For the RN spacetime, we can select the tetrad e^ L as 



Using the ansatz [Toj . 



2M Q 2 \ 1/2 ( 2M Q 2 Y 1/2 \ 

e; = diag ,, 1 _ + ^_\ , h + w >r . rsin£l . (4 ) 



*br,8,4>)=( S^?' (5) 



where for j = I + 1/2, 



l+l/2+m v m-l/2 
(+) _ / V 21+1 Y l 
3 m \ I l+l/2-m v m+l/2 
21+1 Y l 
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and for j = I — 1/2, 



we can get the radial part of the wave equation 



and 



j+1/2— m y-m— 1/2 
21 + 1 Y l 




21+1 Y l 



dr 2 



Vi \F = lu 2 F, (6) 



where 



V\ 2 can be expressed as 



with 



.— + V 2 )G = ^G, (7) 

dr* _ ( 2M Q 2 
dr \ r r 2 

V 13 (r,k)=±^ + W' 2 , (8) 
\ r r z J r r 

where k equals / in the (+) case and — (I + 1) in the (— ) case, and the momentum quantum number I is a non-negative 
integer. One may notice that V\ and V% are supersymmetric partners of the superpotential VF.H'j, and they will give 
the same quasinormal frequencies. So, in the following work, we only discuss the V\ (+) case where the momentum 
quantum number I equals k. 

Fig. ^ shows the potential barrier of the three cases mentioned above with the maximal values 

KnaxG = 5, e = 0.1, Q = l) = 1.441, 
KnaxG = 5, e = -0.1, Q = 1) = 1.689, 

K nax (/ = 5,e = 0,Q = l) = 1.545, 

which can be derived from Eq. (JSJ). And from Fig. we can see the peak of the potential barrier increases with I 
and its position r approaches a limit. 

To evaluate the quasinormal frequencies, we adopt the third-order WKB approximation, which was originally 
proposed by Schutz, Iyer and Will 11] and recently developed by Konoplya to the sixth-order beyond the eikonal 
approximation [13j . This method has wide applications in many black hole cases, and has high accuracy for the low- 
lying modes with n < /, where n and / are the mode number and the momentum quantum number respectively. From 
Ref. , the expression of the quasinormal frequencies is 



\ -,- ( - 2^') 1/2 A] - i I n + \ ) (-2K') 1/2 (1 +5 >i 



where A and ft are the second and third order WKB correction terms 
A = 





" = -77^{ — l$r) (~ — J ) — I -T-77T— ) (' Vl ! "><"" 2 

"(67 + 68^) + ^ (19 + 28a*) 

4a 2 )), (10) 
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with a 2 = (n + 1/2) 2 . Here, Vq is the maximal value of the potential, and V ", Vq", are the second, third and ith 
order derivatives of Vq, with respect to r*. 

Subsituting the effective potential of Eq. (JSJ) into the above formula, we can get the quasinormal frequencies for the 
three classes: e = 0, e < 0, e > 0. To get numerical results, the parameters Q, M, and e should be fixed. Without 
loss of generality, we take M = 1, Q = 0, 0.3, 0.6, 0.9, 0.99, 1 and e = 0, ± 0.1, ±0.3. The results are listed in the 
Tables ffl M Ell and E0 respectively. 

III. DISCUSSIONS AND CONCLUSIONS 

As the main concern of this paper, we will discuss the influence of the charge of the RN black hole on the behavior of 
quasinormal modes in two aspects: how do the quasinormal modes vary with respect to e (the charge of the particle) 
and Q (the charge of the black hole). And following that we will also discuss the momentum quantum number I and 
the mode number n. Thus, this section can be separated into four parts. 

From Eq. 0, we know that we can just consider the sign of the product eQ without taking into account the sign 
of each quantity respectively. As is mentioned above, we take Q to be always positive, and only switch the sign before 
e. Thus, before we move into the discussions, we have to emphasize that positive values of e actually means that the 
product eQ > 0, and vice versa. 

In the following discussions of the quasinormal modes, we focus on the real part and the imaginary part respectively. 
The real part gives us information of the energy of the emitted particles and from the imaginary part we can find out 
the damping time. 

(1) The influence of the particle charge e: From Tabled we can see how the quasinormal modes behave when the 
field changes from neutral to charged cases. 

(1.1) First we study the real part, Re(w): As an example, we study the case, Q = 0.3, I = 1, n = 0. Here three 
values 0.1701, 0.1799, 0.1900 correspond to e = 0.1, e = 0, e = —0.1, respectively. We can see that Re(w) e >a < 
Re(w) e= o < Re(w) e< o. Actually, this result is valid for all cases in TablcQ] In Table ITT1 we can further find out the 
influence of the absolute value of e, and a monotonic behavior is revealed in Fig. [31 i.e., for each given value of Q, 
Re (co) increases as e decreases. 

(1.2) Then we turn to the imaginary part, Im (lu), of the quasinormal modes. Similarly, first we consider the sign of 
e. In Table [IJ we can see that when e is 0.1, the absolute values of the imaginary part of quasinormal frequencies are 
bigger than those of the neutrino field and when e is —0.1, the situation is to the contrary. If we take into consideration 
the value of e, as is shown in Table ITT1 and Fig. we also observe a monotonic behavior, i.e., the absolute value of 
Im (lu) increases as e increases. 

If we bear in mind that the sign of e actually denotes that of the product eQ, we can summarize the above 
conclusions as the following: with larger values of e, the particles have smaller energy and they decay faster. 

(2) The influence of the black hole charge Q: In recent studies reported in Ref. 0, the author investigated the 
influence of the black hole charge on the complex scalar field perturbation (charged) in the RN spacetime. It is shown 
clearly in the figures that the real part of the quasinormal modes of both neutral and charged scalar field grows as 
Q increases, and the latter grows more rapidly. On the other hand, the imaginary part of a given "charged mode" 
approaches the neutral one in the extremal limit (Q = 1). In this paper, following the suggestion of Ref. |8(, we 
calculate the higher multipole number perturbations of the Dirac field with I = 3, 4, 5 and n — 0, and the results can 
be found in Table HJ 

(2.1) The real part, Re (lu): From Fig. [3]we can see it also increases with Q and another character: compared with 
the neutrino field, Re (lu) of the perturbations charged 0.1 and 0.3 are smaller in both the value and the increasing 
speed with respect to Q, however things are different in the cases e = —0.1 and e = —0.3 which shows the influence 
from the Q mentioned above: making the perturbations weaker in the case e > and stronger in the case e < 0. 

(2.2) The imaginary part, Im(w): From Fig. 0]and Table ITT1 we find although no coincidence of the neutrino and 
the charged Dirac field perturbations in the extremal limit exists, there is a turning point separating a monotonically 
increasing region and a decreasing region (note that there is a similar point for the scalar fields, see Fig. 2 of Ref. |£|). 
The explanation of this fact is still unclear but some further work may be helpful, such as finding the relationship 
between the position of the turning point and the value of Q, for different fields. 

(3) The influence of I: It can be revealed by tuning the value of / with the other parameters fixed. As an example, 
taking n = 0, Q = 0.3, e = 0, ±0.1, from Table ITTT1 and Fig. we can see that Re (w) keeps growing with / at the 
speed of 0.2, which shows the contribution of the angular momentum to the energy of the emitted particles. As is 
already discussed above, from Fig. |SJ we can also observe the almost constant difference induced by the values of 
e. From Fig. we find that Im(w) drops sharply as soon as I increases from zero and then it approaches 0.0967 
asymptotically, a value independent of e. This is not difficult to explain. When the particle's angular momentum 
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is big enough, it will hide the influence brought by the different charges. It is also verified by the potential barrier: 
when Q is 0.3 and I is large enough, the maximal potential value with different values of e is always 0.03819^ 2 . 

(4) The influence of n: A related discussion can be found in Ref. [Hj- It can be revealed by tuning the value of n 
with the other parameters fixed. As an example, taking I — 10, Q = 0.3, e = 0, ±0.1, we can see from Fig. [7JIB1 and 
Table HVl that when n increases, Re {uS) decreases and |Im (oj)\ increases. This clearly shows that the low-lying modes 
are of longer relaxation time in not only the schwarzschild spacetime |l0j but also in the RN spacetime, and thus are 
more important for the description of the held evolution around black holes. 

IV. SUMMARY AND SUGGESTIONS 

To summarize, in this paper the quasinormal modes of the spin 1/2 held are discussed at the presence of the 
Reissner-Nordstrom black hole. We study both the real part and the imaginary part of the quasinormal frequencies 
with different values of the particle charge e, the black hole charge Q, the momentum quantum number I, and the 
mode number n. A series of conclusions have been reached, among which the most important one is that, at the stage 
of quasinormal ringing, when the black hole and the field have the same kind of charge (eQ > 0) , the quasinormal 
modes of the massless charged Dirac field decay faster than those of the neutral ones, and when eQ < 0, the massless 
charged Dirac field decay slower. 

Concerning the same problem, the decay of the spin 1/2 field around black holes, there are still a number of 
interesting works to be completed. One work is to consider the massive Dirac field and focus on the effect of mass by 
comparing it with the massless field studied here. The calculation would be more complicated because the mass of 
the particle will appear in the Dirac equation and the effective potential. But we believe that, following the method 
in Ref. [lj}, where the quasinormal frequencies of the (massless and massive) Dirac field perturbations with the 
Schwarzschild black hole are calculated, the massive case with charged black holes can also be worked out. 

One may also be interested in comparing the present results in the RN spacetime with those in some other charged 
black holes such as the RNAdS and dilaton black holes. An especially intriguing case is the Kerr black hole, with non 
zero angular momentum. By comparing the quasinormal frequencies at the presence of non-rotating black holes and 
rotating black holes, like the Kerr black hole|!j, we may understand what influence the black hole rotation brings to 
the quasinormal modes. 
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FIG. 2: The variation of the potential barrier for the massless Dirac held with r at M = 1, Q = 0.3, e = 0.1. 
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FIG. 4: Absolute value of imaginary part of Quasinormal frequencies of e = 0, ±0.1, ±0.3 with Q at M = 1, I — 3, n — 0. 
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FIG. 6: Absolute value of imaginary part of Quasinormal frequencies of e = 0, ±0.1 with I at M = 1, Q = 0.3, n = 0. 
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FIG. 7: Real part of Quasinormal frequencies of e = 0, ±0.1 with n at M — 1, Q = 0.3, I = 10. 
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FIG. 8: Absolute value of imaginary part of Quasinormal frequencies of e = 0, ±0.1 with n at M = 1, Q = 0.3, I = 10. 
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TABLE I: The quasinirmal frequencies for RN black hole, I = 3, 4, 5, n = 0, e = 0, ±0.1, ±0.3, M = 1, Q = 0, 0.3, 0.6, 0.9, 
0.99. 



1 


n 


e 


Q 


= 0.3 




Q 


= 0.6 




Q 


= 0.9 


Q 


= 0.99 








Re(w) 


Im(w) 




Re(w) 


Jm(w) 




Re(w) 


Im{w) 


Re(w) 


Im{w) 


1 





0.1 


0.1701 


-0.1016 




0.1712 


-0.1034 




0.1835 


-0.1018 


0.1879 


-0.0982 









0.1799 


-0.1005 




0.1919 


-0.1012 




0.2207 


-0.0979 


0.2336 


-0.0902 






-0.1 


0.1900 


-0.0993 




0.2138 


-0.0987 




0.2610 


-0.0929 


0.2849 


-0.0808 


2 





0.1 


0.3746 


-0.0974 




0.3844 


-0.0990 




0.4197 


-0.0984 


0.4418 


-0.0936 









0.3847 


-0.0970 




0.4059 


-0.0981 




0.4581 


-0.0963 


0.4891 


-0.0898 






-0.1 


0.3949 


-0.0966 




0.4279 


-0.0972 




0.4978 


-0.0941 


0.5387 


-0.0856 




1 


0.1 


0.3499 


-0.3016 




0.3614 


-0.3059 




0.4008 


-0.3017 


0.4205 


-0.2863 









0.3602 


-0.2999 




0.3834 


-0.3023 




0.4396 


-0.2944 


0.4668 


-0.2740 






-0.1 


0.3706 


-0.2982 




0.4057 


-0.2987 




0.4795 


-0.2868 


0.5152 


-0.2609 


3 





0.1 


0.5726 


-0.0971 




0.5923 


-0.0985 




0.6528 


-0.0976 


0.6912 


-0.0924 






o 


0.5827 


-0 0968 




0.6140 


-0.0979 




0.6915 


-0 0963 


0.7390 


-0 0899 






-0.1 


0.5929 


-0.9651 




0.6360 


-0.0973 




0.7311 


-0.0948 


0.7882 


-0.0872 




1 


0.1 


0.5552 


-0.2953 




0.5762 


-0.2994 




0.6397 


-0.2959 


0.6716 


-0.2799 









0.5655 


-0.2943 




0.5982 


-0.2973 




0.6787 


-0.2913 


0.7242 


-0.2715 






-0.1 


0.5759 


-0.2933 




0.6204 


-0.2952 




0.7185 


-0.2865 


0.7730 


-0.2632 




2 


0.1 


0.5268 


-0.5010 




0.5498 


-0.5076 




0.6178 


-0.4997 


0.6454 


-0.4802 









0.5372 


-0.4992 




0.5719 


-0.5037 




0.6566 


-0.4914 


0.6971 


-0.4578 






-0.1 


0.5476 


-0.4974 




0.5942 


-0.4997 




0.6962 


-0.4830 


0.7445 


-0.4436 


TABLE II: The quasinirmal frequencies for RN black hole, £ 


= 3, 4, 5, n 




0, e = 0, ±0.1, ±0.3, M 


= 1, Q = 0, 0.3, 0.6, 0.9. 


I 


e 




g = o 






Q = 0.3 




Q = 


0.6 


Q 


= 0.9 








Re(w) 


lm(w) 


Re(w) 


Im(w) 




Re(w) 


Im{w) 


Re(w) 


Jm(w) 


3 


0.3 











5525 


-0.0976 




0.5498 


-0.0996 


0.5783 


-0.1002 




0.1 











.5726 


-0.0971 




0.5923 


-0.0985 


0.6528 


-0.0976 









0.5737 


0.0963 





.5827 


-0.0968 




0.6140 


-0.0979 


0.6915 


-0.0963 




-0.1 











.5929 


-0.0965 




0.6360 


-0.0973 


0.7311 


-0.0948 




-0.3 











.6135 


-0.0960 




0.6808 


-0.0959 


0.8130 


-0.0916 


4 


0.3 








0. 


.7488 


-0.0976 




0.7561 


-0.0992 


0.8095 


-0.0994 




0.1 








0. 


.7690 


-0.0970 




0.7991 


-0.0983 


0.8850 


-0.0974 









0.7672 


0.0963 





.7792 


-0.0967 




0.8208 


-0.0979 


0.9239 


-0.0963 




-0.1 











.7894 


-0.0965 




0.8428 


-0.0974 


0.9634 


-0.0952 




-0.3 











.8099 


-0.0961 




0.8873 


-0.0964 


1.0444 


-0.0928 


5 


0.3 











9448 


-0.0972 




0.9623 


-0.0990 


1.0408 


-0.0988 




0.1 











.9650 


-0.0969 




1.0053 


-0.0982 


1.1169 


-0.0972 









0.9602 


0.0963 





.9752 


-0.0967 




1.0272 


-0.0979 


1.1559 


-0.0963 




-0.1 











.9854 


-0.0966 




1.0491 


-0.0975 


1.1953 


-0.0954 




-0.3 








1. 


,0059 


-0.0962 




1.0936 


-0.0967 


1.2759 


-0.0935 




TABLE III: The quasinirmal frequencies 


; for RN black hole, I = 


1, 


5, 10, 20, n = 0, Q = 


0.3, e = 0, ±0.1. 






e 


I = 1 


I = 


2 


I = 


: 3 


I -■ 


= 4 


I = 5 


1 = 10 


I = 20 


i?e(w) 







0.1799 


0.3847 


0.5827 


0.7792 


0.9752 


1.9534 


3.9082 






0.1 


0.1701 


0.3746 


0.5726 


0.7690 


0.9650 


1.9432 


3.8980 






-0.1 


0.1900 


0.3949 


0.5929 


0.7894 


0.9854 


1.9636 


3.9184 


Im(w) 







-0.1005 


-0.0970 


-0.0968 


-0.0967 


-0.0967 


-0.0967 


-0.0967 






0.1 


-0.1016 


-0.0974 


-0.0971 


-0.0970 


-0.0969 


-0.0968 


-0.0967 






-0.1 


-0.0993 


-0.0966 


-0.9651 


-0.0965 


-0.0966 


-0.0966 


-0.0967 
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TABLE IV: The quasinirmal frequencies for RN black hole, n = 0, 7,1 = 8, Q = 0.3, e = 0, ±0.1. 





g 


n = 


n = 1 


n = 2 


n = 3 


n = 4 


jl — 5 


Tl — 6 


n = 7 


Re(w) 





1.9534 


1.9481 


1.9379 


1.9230 


1.9040 


1.8815 


1.8558 


1.8257 


0.1 


1.9432 


1.9379 


1.9276 


1.9127 


1.8937 


1.8711 


1.8454 


1.8171 




-0.1 


1.9636 


1.9584 


1.9481 


1.9333 


1.9144 


1.8918 


1.8662 


1.8379 


Im(w) 





-0.0967 


-0.2904 


-0.4852 


-0.6814 


-0.8794 


-1.0795 


-1.2814 


-1.4852 




0.1 


-0.0968 


-0.2907 


-0.4856 


-0.6820 


-0.8803 


-1.0806 


-1.2827 


-1.4867 




-0.1 


-0.0966 


-0.2902 


-0.4847 


-0.6807 


-0.8786 


-1.0784 


-1.2801 


-1.4837 



